Probing temperature and damping rates in Bose-Einstein condensates using ultraslow light 

experiments 
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We propose a method to probe Landau and Beliaev processes in dilute trapped atomic condensates with a 
multiple state structure using ultraslow light experimental configurations. Under certain conditions, damping 
rates from these collisional processes are directly proportional to the dephasing rates, making it possible to 
determine damping rates through measurement of the dephasing. In the ultraslow light systems we consider, 
Landau decay rates are enhanced at low momenta, which allows one to distinguish between Landau-dominated 
and Beliaev-dominated regimes at the same temperature. Furthermore, the enhancement of Landau rates poten- 
tially provides a way to measure low temperatures (T <C T c ) in dilute condensates more accurately than current 
methods permit. 
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PACS numbers: 

Agreement between mean field approximations and the 
observed behavior of collective excitations in dilute Bose- 
Einstein condensates has been widely verified in low temper- 
ature experiments yj]. However, mean field approximations 
cannot describe many features of higher temperature regimes 
in which damping effects of the excitations become signifi- 
cant. Studying the damping rates of excitations in dilute Bose- 
Einstein condensates allows one to explore thermal field the- 
ories beyond the confines of mean field approximation. 

It has been shown that the damping rates of the collec- 
tive excitations exhibit a strong temperature dependence |2|. 
which has been attributed to Landau processes J5| - the col- 
lision of two excitations to produce an excitation of higher 
energy. Experimental evidence also exists of the reverse oc- 
curring, whereby an excitation decays into two lower energy 
excitations through various processes - collectively known as 
Beliaev processes - such as the coupling of scissor modes |4], 
the decay of Kelvin modes [5], and the decay process signaled 
by the suppression of quasiparticle collisions [6]. There has 
been significant experimental and theoretical interest in these 
damping rates (see for example ||7l lla La 1 1 Of! and references 
therein), yet a fully consistent understanding of these effects 
remains to be achieved. 

In this letter, we propose a method of studying the Lan- 
dau and Beliaev damping rates in dilute trapped atomic con- 
densates with a multiple state structure using ultraslow light 
experimental configurations. Landau processes occurring in 
the multiple-state internal atomic structure of these systems 
are enhanced at low excitation momenta compared to Lan- 
dau processes in systems of atoms with a single internal state. 
This enhancement makes apparent the relative contribution of 
the Landau and Beliaev processes to the overall decay rate as 
a function of excitation momenta, allowing three regimes to 
be identified: 1) dominance of Landau decay, 2) dominance 
of Beliaev decay, and 3) comparable contribution from both 
Landau and Beliaev decay processes. The enhancement of 
Landau dissipation also provides a more effective way to mea- 
sure very low temperatures (T -C T c , where T c is the critical 
temperature) for which, due to the lack of thermal atoms, the 
usual measure of the velocity distribution of an expanded con- 
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FIG. 1 : We consider 3-level atoms illuminated by a probe pulse and 
a classical coupling laser. We assume the condensate atoms to be 
in | bo)- The lasers excite the atoms, into the translational state |cfc). 
7e(fc) represents the decay rate from |cfe). 



densate is not reliable. 

We consider a dilute homogeneous gas of bosons of mass m 
confined in a box of volume V (the effect of the trapping po- 
tential on the homogeneity approximation will be discussed 
briefly below). In the ultraslow light systems we consider, 
three relevant internal states exist: \a), \b) and |c). Where it is 
pertinent to the discussion, we use subscript notation to denote 
the momentum of an excitation within the associated internal 
state. As seen in Figure the ground state \bo), which is a 
weakly interacting Bose-Einstein condensate with No atoms 
in the condensate, is resonantly coupled to excited state \a) by 
a weak probe pulse; and the metastable translational state |cfc) 
of atoms with momentum k is an energy above |6 ) and is 
resonantly coupled to \a) by a classical laser field. Whereas 
in traditional EIT experiments the dark state is a superposi- 
tion of the two lower atomic states and the light fields act as 
classical sources, in the present set-up the dark state is a su- 
perposition of the photons in the probe pulse and the atoms 
in |cfc) fj^foi . and it is the coupling pulse and the conden- 
sate that act as classical sources. The present configuration 
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has been shown experimentally to permit a dramatic slowing 
of light propagation to speeds on the order of 10 m/s 1 14J and 
the storage of light pulses for up to 1 ms before they are re- 
generated 1 15]. (For a recent review see 1 16].) 

If one assumes that 1) the pulse duration is sufficiently long 
that its bandwidth fits within the finite frequency width of the 
transparency window and 2) any incidental couplings to states 
outside this system can be ignored, the dephasing rates, which 
can be measured by the decay of the transmitted probepulse 
energy, will arise primarily from atomic collisions II 2l Il6ll . 
One can show that under typical experimental conditions the 
dephasing rates are proportional to 7 c (fc) fl3ifl7il - the decay 
of the population of the state |cfc)- by a proportionality factor 
close to 1.0 (slightly greater for light delay experiments than 
for light storage experiments) [13]. 

The population of state \ck) decays mainly due to 2-body 
elastic spin-conserving collisions with the atoms in internal 
state |&) characterized by a positive scattering length at, c . To 
leading order, collisions between atoms in \b), characterized 
by a positive scattering length abb, change the free particle 
dispersion relation to the Bogoliubov quasiparticle dispersion 
relation, e h p — [ip\/2 + p 2 111 (811 . where [i = (hko) 2 /2m is 
the chemical potential and p is the normalized excitation mo- 
mentum. (Unless otherwise specified all momenta are nor- 
malized by the healing length fc^ 1 = l/y/Snnoabc where 
no = No/V.) Higher order collisions among atoms in \b), 
which cause the quasiparticle populations to decay, can be ig- 
nored in the dark state picture because the atoms in \b) serve 
only as a classical source in our system. To allow us to treat 
the atoms in \b) as if they are in thermal equilibrium, we as- 
sume that atomic collisions within internal state \b) occur at 
much shorter time scales than atomic collisions between \b) 
and |cfc), i.e. abb S> ab c 1131 . Collisions among atoms in in- 
ternal state |c), characterized by a positive scattering length 
a cc , can also be ignored as there is no macroscopically oc- 
cupied mode in |c). Consequently the dispersion relation for 
excitations in \c) is the free particle dispersion relation given 
by £p = V-V 1 (note that the total energy of excitations in |c) 
is e"p + ebc) 12111 . We also assume that abba cc > a\ c in order 
to avoid phase separation 1 22] . Finally, collisions involving 
atoms in \a) can be ignored because the population of \a) is 
negligible. 

Assuming the system of atoms is near thermal equilibrium, 
one can determine 7 f1 (AQ by c alculating the imaginary part of 
the perturbed energy I23H24I1 . by using the Golden Rule fioll 
or by taking a time-dependent approach 013. The decay 
rate of \ck) comprises the Beliaev and Landau decay rates, 
i.e. 7c(&) = lf(k) + 7^ (k). The Beliaev decay rate, cor- 
responding to the collision between a particle in \ck) and a 
condensate atom in |6o) to produce one quasiparticle in \bj) 
of energy e b and one particle excitation in \ci) of energy e£ 
such that ei = ef 



A, is given by 



where the interaction strength for the contact potential be- 
tween atoms in \b) and |c) is given by Ubc — 4:irh 2 ab c /rn 
where n b - = — . b .] . and n c „ = — rr — - — 1 tw , „, 1 

are the thermal populations associated with \bj) and |cj), re- 
spectively. The Landau decay rate is caused by the reverse 
process - a particle in |cfc) collides with a quasiparticle in 
\bj) to produce a particle excitation in |cj) with higher energy 
e i ~ € k + e 5 anc ^ a condensate atom in |6q) - and is given by 
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where nl — 



{ui-vtfM-nZMG + e J-ef) (2) 



(Note that the differ- 



exp[(e|+e bc +ep/fc B T]- 



ence in n c i between the two cases will become irrelevant, as 
discussed below.) Uj and Vj are the quasiparticle amplitudes 
that characterize the Bogoliubov transformation 1 18] and the 
square of their difference is given by (u p — v p ) 2 = p 2 7 ' z h p 
where the dimensionless energy scaled by the chemical po- 
tential is given by z h p — +p 2 . 

As we are considering a homogeneous system, the sums can 
be converted into integrals to arrive, after some manipulation, 
at 



left) = yfno 




1 



(l+n?+nf), 



(3) 



and 



dz° 1 



(nS-nf). 



(4) 

Since ebc/ksT ^ 1 in typical experimental situations in- 
volving condensates (note that we assume ebc/ksT ^> 1 in 
this paper unless otherwise specified), the number of ther- 
mally excited atoms in state |c) is negligible, i.e. n\ « 0, 
implying 



7 c L (fc) 



n » 3 M/(r) 



(5) 



(1) 



where the dimensionless temperature is given by r = fcsT//z 
and J(t) is a function only of temperature and is given by 

J °° dz b ^1 — — j==g==\ nj(r). This assumption of a large 

energy gap between states \b) and |c) implies a large differ- 
ence in the thermal occupation numbers in \b) and |c), which 
has a significant effect on the Landau decay rate of the low 
momentum excitations in |c), which now diverges as 1/k as 
k — > (the assumption that the excitations have well-defined 
momenta breaks down at small momenta, an effect that we 
describe below). Compare this to the single state case where 
7 L (fc) — * as k — > 0, e.g. 7 L (fc) oc fcr 4 for k < t < 1 and 
7 L (fc) oc kr for t > 1 and k <C 1 (see J26ll and references 
therein). At high momenta where k 1, the Landau decay 
rate falls off as 1 jk as in the single state case l27il . 
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FIG. 2: The inverse of the decay rate as a function of excita- 
tion momentum k (scaled with the healing length) for T = 0.1 T c 
(r = 0.86) for the upper plot and T = 0.2 T c (r = 1.72) for the 
lower plot. The solid line represents l/ n / c (k) - the inverse of the to- 
tal decay from \ck) assuming the experimental realistic scenario of 
£i)c/fcsT ;§> 1. The dashed line represents the inverse of the decay 
due to Landau processes, l/y^(k), and the dot-dashed line repre- 
sents the inverse of the decay due to Beliaev processes, l/j^(k), 
both assuming et c /fcsT ;§> 1. To illustrate the importance of etc 
we include l/~/ c (k) as a dotted line assuming etc = 0. The gray 
area signifies the excitation momenta that can be measured before the 
edge of the trap is reached, i.e. l/^ c {k) < h 2 fc °^ . We assume typi- 
cal experimental parameters: a = 3 nm, T c = 500 nK, fj,/h = 2-7T 
kHz, ko — 1 /im" 1 , R = 100 fim (along the major axis of elongated 
condensates). 



The enhancement of Landau decay rates at low momenta 
in this system allows one to determine whether the Landau 
or Beliaev decay is dominant: the Landau decay rate de- 
creases monotonically with momentum (unlike in the single 
state case) whereas the Beliaev decay rate increases monoton- 
ically with momentum as seen in Figure I3 This allows one 
the possibility to independently probe Beliaev- and Landau- 
dominated processes at the same temperature. Eventually, as 
the temperature increases, the Landau processes will become 
dominant over a wider range of momenta because Landau pro- 
cesses can involve any mode with an appreciable population 
whereas Beliaev processes are restricted to modes whose en- 
ergies are less than that of decaying particles. 

We have assumed a homogeneous gas, which should be 
a good approximation for the central region of the large 



alkali-metal condensates achieved in present-day experi- 
ments. However, we have implicitly assumed that excited 
atoms in state \ck) remain in the interaction region, which is 
ultimately limited by the size of the condensate (for large con- 
densates this is the Thomas-Fermi radius), for times longer 
than that of the characteristic decay, i.e. 7 c (fc) > 2fc-^-^-^ 
where R is the characteristic size of the condensate. This puts 
an upper limit on the excitation momentum one can measure; 
the measurable momenta are represented by the gray region in 
Figure|21 

We have implicitly assumed that the excitations have well- 
defined momenta, i.e. j c (k)/el <C 1. For the Landau- 
dominated regimes at low momenta this condition implies that 
k 3 v / 87TV / ^o^fc/( T )' which is not very restrictive be- 
cause we assume a dilute gas such that y/ «o fl b c <C 1> an d 
/(t) < 1 for temperatures in the range of r < 1 where both 
Beliaev and Landau processes are relevant. For the Beliaev- 
dominated regime at relatively high momenta (k > 1) and 
temperatures in the range of r < 1, the zero-temperature re- 
sult -ff(k)x=o = g(k)\/noal c n/h dominates; here, g(k) = 
%/87rfc[l - ln(fc 2 + vT+ k 4 )/k 2 ]. Th e assum ption of well- 
defined momenta is satisfied when \/noa^.g(k)/k 2 <C 1, 
which is always true for a dilute gas since g(k)/k 2 < 1 for 
all k. 

The enhancement of the Landau decay rate at low momenta 
has another important consequence: it offers a method of ac- 
curately determining very low temperatures, i.e. T -C T c . At 
T = the Beliaev damping is the only remaining dissipa- 
tive process, as is the case with the single-level system, and 
7c(fc)T=o = jf(k)T=o — ► as k — > as k 5 (the full expres- 
sion is given above). Unlike in the single atomic state case, 
however, at low temperatures r = fc^T/ /i <C 1 Landau decay 
dominates at low momenta and can be written as 

7cW-7 c L W-^^V^L^ (6) 

where ( is the Riemann zeta function and the excitation mo- 
menta remain well defined for fc 3 3> y / riJ)a|^T 3 . 

The decay approximation used in this paper is only strictly 
valid for our homogeneous system with its continuous spec- 
trum of excitations (and perhaps also for anisotropic parabolic 
traps where the discrete energy levels are chaotically dis- 
tributed 1 9]). In experimental situations using a trap, however, 
the spectrum of excitations becomes significantly discretized 
at lower momenta, which would considerably impede Beli- 
aev processes at low momenta. Nevertheless the argument 
remains unaffected that Landau decay rates are enhanced at 
low momenta due to the large energy gap between \b) and |c). 
We maintain in this paper that this 1) makes it possible to dis- 
tinguish between the Beliaev and Landau decay rates at the 
same temperature and 2) provides a sensitive low temperature 
thermometer. 

The author gratefully acknowledges stimulating discus- 
sions with David Guery-Odelin and Sam Morgan. 
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